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Q ! Abstract 

OS . 

. In a previous paper || , we derived a post-Newtonian approximation to cosmology which, 

in contrast to former Newtonian and post-Newtonian cosmological theories, has a well-posed 
Oh initial value problem. In this paper, this new post-Newtonian theory is compared with the 

fully general relativistic theory, in the context of the k = Friedmann Robertson Walker 
cosmologies. It is found that the post-Newtonian theory reproduces the results of its general 



X 



OA | relativistic counterpart, whilst the Newtonian theory does not. 

1 Introduction 

0^ , The equations of general relativity are difficult to solve and often become tractable only in a 
Newtonian context. Hence, it is desirable to use Newtonian theory rather than general relativ- 
ity where possible ||. By Newtonian, we mean that theory which is obtainable from general 
relativity by taking the limit of weak gravitational fields and small velocities. Unfortunately, 
when Newtonian theory is applied to cosmology, the boundary conditions become redundant, 
and subsequently the Poisson equation no longer has a unique solution. Thus, the theory is no 
longer well-posed and causality may be violated. 
(2jT)! In a previous work Q , we explored the insufficiencies of the Newtonian theory in detail and 

showed that in addition to the lack of well-posedness Newtonian theory is also incomplete. The 
reasons for these insufficiencies are the following. In Newtonian theory, which is the expansion 
of general relativity in terms of weak gravitational fields and small velocities up to order c -2 , all 
that remains of the field equations is the Poisson equation, ie. the evolution equation and the 
constraint equations have become one and the same. With no time evolution equation for the 
Newtonian potential eft, a unique solution for <p cannot be found, and so the theory is not well- 
posed. Furthermore, the term - y^), a term of order c -4 , does not feature in the field equations, 
and hence the Bianchi identities cannot be fully obtained. Although completeness is obtainable 
at order c -4 , it is necessary to go to order c~ 6 and then to reformulate the field equations 
as wavelike equations, in order to obtain a well-posed initial value problem. Throughout the 
remainder of this paper we will call this approximation to general relativity "post-Newtonian 
theory" . 

In the following we will explore the Friedmann Robertson Walker (FRW) cosmology in the 
context of the above Newtonian and post-Newtonian approximations. It has been shown that 
the FRW models have a simple Newtonian interpretation |7|, GJ. Howsoever, it will be outlined 
here that the Newtonian theory cannot reproduce all the solutions of the fully general relativistic 
theory: Variations in the equation of state do not cause changes in the solutions of the theory 
since the pressure does not enter into the dynamics. This is reflected in the Raychaudhuri 
equation being reproducible only for the case of vanishing pressure. Thus, the Newtonian theory 
is only useful for the special case of dust. In the post-Newtonian theory, on the other hand, such 
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difficulties are overcome. The pressure does feature in the dynamics, and varying the equation 
of state produces correspondingly varying solutions, allowing for the full range of possibilities of 
its general relativistic counterpart. We therefore argue that the post-Newtonian theory should 
be used whenever we have non-zero pressure in the univese. 

In Section 2 we write down what constitutes the general relativistic k = FRW cosmology. 
In Section || we will introduce the Newtonian and post-Newtonian approximations and put the 
FRW metric into a form such that we may draw comparisons between these three theories. In 
Section ^| we explore the Newtonian approximation, and we will show that the theory may only 
be used in the special case of dust. In Section || we will find that the post-Newtonian theory is 
able to fully reproduce the results of the general relativistic case. We will discuss our results in 
Section |6[ 

2 The Fully General Relativistic k = Friedmann Robertson 
Walker Cosmology 

In the case of the flat FRW metric, the field equations are the Friedmann equation and the 
Raychaudhuri equation 




(1) 



and 




4ttG(p + 3 P c- 2 ), 



(2) 



from which one may obtain the Bianchi identity 




= M,R~^ 1+W \ with C a 
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Typical values for w include: 



Matter (w 



0): 



p oc R 



-3 



R oc t3 , 



Radiation (w = — 



p oc R 



-l 



R oc ta, 



Stiff Matter (w = 1) : 



p oc R 



6 



R oc ta. 



(3) 
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3 Newtonian and Post-Newtonian Approximations and the 
Friedmann Robertson Walker Metric 



We now would like to study the FRW cosmology in the Newtonian and post-Newtonian approx- 



imations. Following a scheme similar to that of Weinberg's [11], we adopt units in which the 
typical velocity has magnitude 1, i.e. ~ c" 1 , and assume a one parameter family of metrics 
9nu{x x ,c) for which there is a system of coordinates (x°,x l ) in which the components of the 
metric have the following asymptotic behaviour as c — ► oo: 



5oo 

9ij 



-1 - 20c -2 - 2ac -4 - 2a'c -6 



2a" a' 8 . 



c 



2(j)5ijC 2 + ctijC 4 + ol iA c 6 + a''c 8 



(4) 



The usual Newtonian theory is obtained as the 0(cT ) limit of (Q), while the complete 
Newtonian approximation is the 0(c~ 4 ) limit. Reformulating the field equations of the 0(c~ 6 ) 
limit as wavelike equations defines the post-Newtonian theory (see ||]). 

We need to put the k = FRW metric, 



ds 2 = —dx\ + Ej-R(i) 2 dxf , where xq = ct, 

into a form of which we can read off the potentials (f>, d, a and aijj. To do so, we consider the 
following coordinate transformation 



x 

X i 



Tc + tc' 1 +r'c" 3 , 
R~ x Xi + X iC~ 2 + XiC-\ 



with 



T 
/ 



Xi 

x\ 



A^ + Aij^Xij, 
Bij(t)Xij + Bijki(t)Xijki, 
Cij(t)Xj + Cijki{t)Xjki, 

Dijkl{t^)Xjfci + Dijfci mn (t)X jklmni 



where Xij = X^Xj and similar for X^ etc. Throughout the remainder of the paper we assume: 
A <C c 2 to ensure convergence of our expansion in c~ 2 . A(t), Aij(t), Bij(t), Bijki(t), Cij(t), 
Cijki(t), Dijki(t) and Dijki mn {t) are arbitrary functions of time. In these coordinates the metric 
becomes 



ds' 



c 2 dT 2 



-l + c~ 2 {-2A-2A ij X ij + 



R 
R 



Xj. 



-c 4 (—2BijXij — 2BijkiXijki + (A + AijXij) 2 — 2RCijXij — 2RC i j] t iX i j] t i 



+cdTdX l 



4A ij X j -2^X,, 
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-c 3 (—4BijXj — SBijkiXjki — 4 yj^j AijXkkXj + 2RCijXj + 2RC; L j k iXj k i 

R R 

-2RCijXj — GRCijkiXjM + 2—AXi + 2—Aj k Xij k ) + 0(c 



+dXidXj 



R 

5ij + c~ 2 (-4A ik AjiX k i + 2i?Cj, + 6RCijkiX k i + 4—A ik Xj k 

K 



+c 4 (—8AA ik AjiX k i — 8B ik AjiX k i + 8A ik AjiA mn X k i mn — 16AjiB ikmn X k [ mn + 6RDij k iX k i 
+WRDij klmn X k i mn + _R (Cy + 3Cij k iX k i) + 4 I — 1 A ik AjiX k iX mm + ARAuCj k Xi k 

R 

—4RA in Cj mk iX mnk i — ARA a Cj k Xi k + \2RA in C j mk \X mnk i — 4—AAuXji 

R R R 

—4—A im A k iXj k i m + A—B ik Xj k + 8—B imk iX mk ij) + 0(c~ 6 ) 

XI It It 

Since in (0) there are no terms of order c _1 , we require that 



(5) 



(6) 



In addition, the terms of order c 2 in goo and gij may be identified, and from the powers in 
Xj one can read off 



and 



A 
R 



a 



ijkl 



4 The Homogeneous and Isotropic Newtonian Cosmological 
Theory 



The Newtonian theory is the c 2 cut-off of (|j) and consists of the field (Poisson) equation 



AirGp, 



(7) 



((f) is the Newtonian potential and p is the density) and the continuity and Euler equations 
of fluid dynamics 



P + pv-i,i = 0, 

Vi + 4>,i + -p,i = 0, 
P 



(8) 
(9) 



where Vi is the velocity field and p is the pressure. Homogeneity implies that the density and 
pressure are merely functions of time and that the velocity field is the same relative to all 



observers. It can be shown, || and [10], that this amounts to V{ = Vij(t)Xj. By substituting 
the Poisson equation into the Euler equation we see that the Newtonian potential must be of 
the form <f> = aij{t)XiXj + a(t). Therefore, the Newtonian approximation of a homogeneous 
cosmology is 



an = 4nGp, 
P + pVu = 0, 

Vij + V ik V k j = aij. (10) 

The FRW cosmology is the most general isotropic and homogeneous solution. Thus, we will 
only consider the case where ( JlO) ) becomes isotropic, ie. shear-free and rotation-free. Heckmann 
and Schiicking in || formulate the more general anisotropic Newtonian cosmology where there 
is shear and rotation. To this end we make the following decomposition 



Vij = ^9Sij + + Wij, (11) 



where 



= Vu, 

1 

= e jik w k = -(Vj - Vji), 

The trace part 9 is the expansion, the tracefree symmetric piece &ij is the shear and Wij, 
the anti-symmetric part, is the rotation. Using this decomposition in the Euler equation and 
setting a = and w = 0, we get 



±6 2 -4ttG P , 

T'ij = Ti^kk^ij j 



with continuity equation 



p + P 6 = 0. (12) 

Defining a function R'(t) such that 9 = 3^; then the solution of the continuity equation 
(0) is 

p = C'R'- 3 , (13) 

with C a constant. Finally, using ( |l3|) in the Euler equation, we see that the Newtonian, 
isotropic, homogeneous cosmology is given by 



R_ 
R' 
P 



A-rrGp, 

~3' 
C'R'-\ 



-irGp, 



(14) 

(15) 
(16) 
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Notice here the use of R(t) as the general relativistic scale factor in (|l|) through to (||), 
and R'(t) as the Newtonian scale factor in the Newtonian theory ( O ) through to (16). We are 
now in a position to compare these two theories. The general relativistic theory is well-posed. 
Equations (|l|) and (§) are consistent with the Bianchi identities. In the Newtonian theory there 
is only one field equation (|l4|), and there is no completeness because ( |l4| ) does not give ( |l5| ) 



and (16), and nor is the theory well-posed. Also, notice that in the general relativistic theory 
pressure occurs in the dynamics of the theory, whereas in the Newtonian theory, pressure does 
not occur anywhere in the dynamics and is only defined through an equation of state. 



Equation (15) has the same form as the Raychaudhuri equation (at least when p = 0). Using 
the Raychaudhuri equation it may be deduced that R' oc VS. Thus the equations (|l5|) and ( |l6| ) of 
the Newtonian theory predict the same results, at least for the case of matter, that the general 
relativistic equations (||) and (||) do. 

How, if at all, do R(t) and R'(t) differ? To answer this question we use the remaining piece 
of information - (|l4|), the Poisson equation. Since <f> = aij(t)XiXj + a(t) is the term of order c~ 2 
in the c 2 dT 2 piece of the FRW metric (|E|); a comparison of (Q) and (0), and using (|6|) yields 

a ii — R ' 

and 

a{t) = A, 

from which, with the aid of (|lj|), we can deduce 

I = -r ?- < 17 > 

This is again the Raychaudhuri equation of general relativity for the case of vanishing pres- 
sure. (Since the Newtonian potential (j) only appears as 4>^k = 2a&fc(i) in equations (Q) to ([!]) 
we may set the piece a{t) = without loss of generality; then A is just a constant.) 

Hence, the correct general relativistic scale factor R{t) is equivalent to the Newtonian scale 
factor R'(t). Dautcourt [g], shows how the validity of the Friedmann equation within Newtonian 
cosmology can be understood: Newtonian cosmology is applicable only when confined to a 
neighbourhood of the observer, corresponding to distances which are small compared to the 
Hubble distance. 

Although reproducing results similar to the general relativistic theory, the Newtonian theory 
suffers in that varying the equation of state will have no effect on the outcome of the solutions 
for p(t) and R'(t). This being due to the fact that the pressure has not appeared anywhere in 
the dynamics. Thus we can only reproduce the results of the matter dominated case of general 
relativity. 

5 The Homogeneous and Isotropic Post-Newtonian Cosmologi- 
cal Theory 

The field equations for the post-Newtonian theory || are 

(p M = 4nGp+ ^(-fefe-2t), (18) 

d,kk = l6irGpvi + ^ (J)ijj - »j) . (19) 
4>ij - c 2 4> ijtkk = T>y + c 2 [l6TrG(pviVj + <%p) - %] , (20) 
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where the a and a^- of (01) are such that 



kj = a-ij - 25ija, 



with 



21 



23; 



60,^,1 - 167rG(,9t; 2 + 4p0), 
•X.,,/0.> + 2Cj0,ij - 20jCi,i 



16ttG 



1 



VjP + pViV 



(&(k,kj + CjC 



$ij(6<f>,k<t>,k + 327rG/?0), 



Oc (Cijfc + Cj.ifc) + 2(k(k,ij + (k,i(k,j + (i,k(j,k 



k,ki J 

-20,* (0M,j + 0jy,i - 20^, fc ) - 1600,j0,j + (p^cpkkj 

+<t ) ,j ( t ) kk,i - 20 {4>ki,jk + 0fcj,ifc — 0ij,fefc — 4>kk,ij) — 24>ki4>,jk 

-24>kj4>,ik + 24>,ij4>kk - Sij 

—(k(m,mk — 4:4>,k(f>km,m + ^<\>,k$mm,k ~ 1200 *,fc0 \k 
^km,mk - 2(/) mm ,kk)} + 8nG [2pViVj 

-(bp 



2 Ck,mCk,Tn ~i~ ^Qm^Ck^m ~\~ ^ (Cfe,fe) 



+2p(20 + U + p(j)ij + 5jj 2p0u 



3 _, 1 j. 

+^9,kP,k + ^P9,kk 



Along with the harmonic gauge conditions 



(21) 
(22) 



this system forms a closed set which is consistent because the Bianchi identities are obtainable 
from the field equations: 



and 



+ (P V j)j ( 1 + „ 2 



IT 



+ 



p ( 2vjVj + 2vjV k v kJ + -Cj,j 



{2'ft,iCi,jj 2( l i(/)jii 3£j,j0j'j) 



P{vi + VijVj + 0,j) + P 

23 



1 



-(2l + S3 i ,>;+ 



16vrGc 2 

QSijj - 2^(21,1 + 4>jk,kij) ~ + fijkjk) 



(23) 



(24) 



Specifically, (^) and the , i derivative of ( |19| ) imply 3/c?i of the first constraint equation 
fll8| ) , while (^4|) and , i of eq. (^0|) imply the time derivative of (|l9|) . 

To make comparisons with the FRW cosmologies we need to consider the case where this 
theory is both homogeneous and isotropic. The following are the most general ansatze, providing 
isotropy and homogeneity, for tensors of rank n, expanded up to order X n+2 : 
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= &<Kt)x 2 + 

Q = W((t)X 2 X l+ V>at)Xi, 
0y = ( 5 )0(t)X%-+ ^{t)X 2 Xij+ ^^(t)X\ j + ®$(t)X ij + (1) 0(i)<%. (25) 

Again, we decompose the velocity field, Vy, (|TTJ) and set Wij = and <Jij = 0. Thus we have 
Vi = VijXj = \0'Xi. Pro tempore we do not make the substitution 6 = 3^-, instead, we put 9'. 



The field equations, fllq) to (pQ), are then 



3 < - 2 U = 2nGp + c- 2 



10 ( 5 >0 + 15 ( 4 )0 + 12 ( V - SirGp Q#' 2 + 4 ^0 1 0. 



5 < 2 >C - *irGp\# + c- 2 -( 3 )0 - 2 ( 2 )0 + 9 ( 2 )0 «C + 8^G ( -~0'p + ~p «C 



1„/ , 1 



0. 



2 (5)£ + 3 (4)^ _ 6 (2)^(2) C + 8 ^ Gp f 1 _ I (2) C \ = 0> 



V27 2 



6 (3) + 2 ( 2 )0 _ 16 (2)0 (i)^ + i6vrG(p + 2/) W0) 

.(i) J + H (i) C 2 - 4 (2)0 (D0 + 4 W0 (2)0 + 8 «0 ( 3 )0 + 8nG L{ «0 + 3 {2 U) ~ \p (1 U 



20 ( 5 >0 +2^0+8 (2)0 2 + 32nGp ^0 + c~ 2 + 11 «C (2) C - 44 ( 2 )0 ( 3 )0 

-16 (2)0 (2)0 + 32 «0 ( 5 )0 - 8 «0 ( 4 )0 + 48 (2)^ (D0 + (p(^4> + 2 (1) 0^' 2 ) - - (2 )# 

V 9 2 



14 W0-16 (2)02 + IGnGp-e' 2 



+c 



.(2) J + 2 W( (2)C - 64 (2)0 2 (1)0 - 8 (2)0 ( 3 )0 - 16 (2)0 ( 2 )0 - 16 (1)0 ( 5 )0 + 4 «0 ( 4 )0 



+8vrG ( 2\e a (p + 2p ^0) + p ( 2 >0 



0. 



(5) J + ^ (2) C 2 + 60 (2)0 (5)0 _ 2Q (2)0 (4)0 _ 4§ (2) 3 _ ^ Gp ( 2 (2)0^/2 + (5) j = ^ 

2 V 9 



(4)0 _ 6 ( 2 )C 2 + 64 ( 2 )0 3 - 32 ( 2 )0 - 12 (2)0 (4)0 
-SvrGp ^4 ^0 V + 2 ^' 4 +(4) 0) = 0. 

The homogeneous harmonic gauge conditions (|2l|) and become 



«0 



4 

.£ (2)A 

4 



Using ( p8|) and the time dervivative of (|26|), the harmonic gauge condition (^) reads 



2 (3) + 4 (2 )0. 

4 (5)0 + 6 (4)0 5 



(35) 

(36) 

(37) 
(38) 



p + e'( P + pc- 2 ) + c 



.^(2)0« C -4p( 1 )0 + ^( 1 )C 



0, 



(39) 



where (27) and ( |29| ) have been used to express the c~ 2 parts of the equation in this relatively 
simple form. This is the continuity equation ( |23"|) for the special case of homogeneity and isotropy. 

Next, consider equations ( |3l"D and (|32|) with the time derivative of (|2"8|). These equations, with 
the help of ( |33|) and (|34j) , provide us with the Euler equation, fl24|), for the case of homogeneity 
and isotropy: 



,(^ + ^ +2 ^) + _l_, 



265 (2 )0 2 «0 + 132 (2 )0 ( 3 )0 + 120 (2 )0 (2)0 



-18 (2)0 «C - 18 (2)0 «C - 30 «C {2) C + 80 (V 5) <^ + 120 ( 4 ^ + 54 ( 2 )0 «C V 

3 

+16.G (-p Q w c + w^) + pV +p + ^ (2 V - 7 V 2 ) 

+p _(3) _ 14 (D0 V _ 2 (2)0 _ 24 (2)0 (1)0 _ 8 (1)0 \ e+ l _ (D^l = 0. 

\ 2 9 3 2 y y 

(40) 

We may now read off the potentials (/>, Q an d from the FRW metric (|5|). With the help 
of the harmonic gauge conditions ( |35| ) to (|38| ) we are able to write down the functions ^4>, 
etc. which appear in the ansatze (p5[): 



«0 

(2)0 
(D C 

(2)^ 

«0 
(2)0 



A, 

1 R 

~2R' 
4 

~3 J 
2 ('R _ RR\ 
b\R R? ) ' 
4i 2 , 

2 ... 8 -R 2 - R ■ 
— A H A— + -A— + 2A 

5 15 R 



R' 
R 
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(41) 

Using these potentials in (40) and adopting a power expansion for 0' 



+ c- 2 6", where = 3$, 



one obtains, with the help of (|30| ) and after considerable simplification 



-4vrG(p + 3pc 



+ c 



;R 



R 



3A 



R 



(42) 



with the required inclusion of the 



This is the Raychadhuri equation of general relativity 
pressure term, which was missing in the Newtonian case 
i < c 2 , we can be sure that the corrections are of higher order and do not contribute to the 
theory at the c° level. We may now use this result, along with the now known values for the 



17]), and with c corrections. Since 



potentials (41), in ( |39| ) to derive the Friedmann equation 



8irG 



P + Te- 



rn 



where T is defined to be the solution of the differential equation 



2 R\ R . 

•5 5 +6 5 u 1 



RR 




+ 20" | i| 

3R 




The Bianchi identities (|39|) and ([4Q), the field equations (|2J|) to ( |34[) and the harmonic gauge 
conditions (|3~5| ) to (38) may be understood in the following manner: The time derivative of 
equation fl28|) along with the equations (31), ([32]), (|33|) , (134]) and with the aid of the harmonic 
gauge conditions give us the Bianchi identity (^) which leads to the Raychadhuri equation 
(|4"2|). The time derivative of equation (^) with equations (|28|), (p7|), (pS|) and the harmonic 
gauge condition (|3^) gives the Bianchi identity (|39[), which leads to the Friedmann equation 
d43|). This leaves us with the field equation ( p0| ) which may also be combined with (^) to give 
the Raychadhuri equation. 

It may also be shown that the following relationships between the field equations exist. The 
time derivative of (^) is (|28|). The time derivative of (27) with the aid of (|30| ) is (|29|) . Equations 
( |3l| ) and p2] ) may be combined in such a way as to be (^). Finally, Equations (|3^) and ( p4| ) 



may be combined in such a way as to be the second time derivative of (£ 

Thus the post-Newtonian approximation may be completely defined for the FRW cosmology 
with the Raychadhuri equation (^), the Friedmann equation (|43| ) and an extended Poisson 
equation (p6|), which after substitution of (|4l|), becomes 
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3^ = -4irGp + c" 2 (- A -lfaGpA). (44) 
R 

The unknown 6" may be set to zero without loss of generality. Equation (^) may be 

equivalently represented by (|3^). Given an equation of state (pc~ 2 = wp, w an arbitrary 



constant), the system (39), (|42|) and (Q) gives equations for p, R and A respectively, and thus, 
forms a well-posed set, allowing the unknowns R(t), p(t) and A(t) to be determined. 

It is the gravitational potential A{t) that incorporates the pressure into the theory. The 
pressure and A{t) are related through a third order differential equation one can obtain by 
combining (^) and (|4"4|). Varying the equation of state therefore results in varying solutions 
for the density p(t) and for R(t), which are analogous to those solutions predicted by the fully 
general relativistic theory for the FRW cosmologies, (||). Setting A = 0, for example, implies 
vanishing pressure, p = 0, which is the post-Newtonian approximation of FRW for dust. In this 
case the higher order corrections of the post-Newtonian theory disapear and one ends up with 
the Newtonian theory, which we have shown in Section [| describes the case of dust. 



6 Concluding Remarks 

The FRW cosmology is a very good approximation to the large scale structure of the universe, at 
least to the present epoch Q. Due to its obvious simplicity, Newtonian approximations, where 
they reproduce results which are similar to the fully general relativistic theory, are preferable. 
General relativity has a well-posed Cauchy problem in the case of perfect fluids with a barotropic 
equation of state ||. Newtonian theory, on the other hand, is not well-posed ||. We have seen 
that the Newtonian theory reproduces the results of it's general relativistic counterpart only for 
the special case of dust, ie. for a matter-dominated universe with an equation of state, p = pc 2 . 
Although, any equation of state may be written down, it won't make any difference since the 
pressure does not appear in the dynamics. 

The post-Newtonian theory, which is consistent and well-posed, does provide field equations 
with pressure entering the dynamics through the potential A(t). We are thus able to vary the 
equation of state and, in doing so, will obtain various solutions for R(t) and the density, p. 
Hence, the post-Newtonian theory seems to be a favourable approximation of the fully general 
relativistic theory. 

I would like to acknowledge Peter Szekeres for useful discussions. 
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